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We address the problem of determining the type I, type II or borderline dual superconductor 
behavior in maximal Abelian gauge SU{2) through the study of the dual Abrikosov vortex. We 
find that significant electric currents in the simulation data call into question the use of the dual 
Ginzburg Landau Higgs model in interpreting the data. Further, two definitions of the penetration 
depth parameter take two different values. The splitting of this parameter into two is intricately 
connected to the existence of electric currents. It is important in our approach that we employ 
definitions of flux and electric and magnetic currents that respect Maxwell equations exactly for 
lattice averages independent of lattice spacings. Applied to specific Wilson loop sizes, our conclusions 
differ from those that use the dual GLH model. 
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I. INTRODUCTION 

In the search for simplicity in the physics of color con- 
finement, lattice gauge theory models have been stud- 
ied extensively for a clue to a mechanism or an under- 
lying principle governing the phenomenon. Spontaneous 
gauge symmetry breaking (SGSB) of the dual £7(1), and 
the resulting condensation of U(\) monopole currents, 
defined after appropriate gauge fixing, remains a candi- 
date. The persistent monopole currents of dual super- 
conductivity in pure U{1) lattice models leads to con- 
finement of charge. There has been some success in the 
postulate of Abelian dominance in correlating monopoles 
and confinement physics but no breakthrough in uncov- 
ering a definitive mechanism. There are a number of 
reviews of the subject [1-9]. Some of the principal di- 
rections include accounting for string tension in Abelian 
Wilson loops [10, 11], similarly for monopole dominance 
of Wilson loops[12, 13], correlating percolating monopole 
clusters and confinement [14-17], spontaneous gauge sym- 
metry breaking (SGSB) of dual L7(l) symmetry of the 
vacuum[18-20], the dual Ginzburg Landau Higgs (GLH) 
model[21, 22] and studies of the structure of the dual 
Abrikosov vortex in the confining string [23-31] modeled 
by the dual GLH model. 

Truncation to relevant variables invariably leads to sys- 
tematic errors. Ambiguities due the requirement of fix- 
ing the gauge [32-36] and Gribov copies [37] contribute 
to these perennial problems. Further it is difficult to see 
how this mechanism alone can explain all aspects of color 
confinement in arbitrary systems. Nevertheless one can 
argue that the well established phenomena can be part 
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of a larger picture. 

One possibility is that we are seeing SGSB as a general 
principle rather than just a confining mechanism. The 
Pisa group[18-20] has given extensive evidence for this 
through QCD vacuum studies of the deconfining transi- 
tion. SGSB also manifests itself in the existence of a dual 
Abrikosov vortex in the confining string[23-31]; a com- 
plementary approach giving rather direct evidence for the 
same principle which is the subject of this paper. 

The dual Ginzburg Landau Higgs model has 
provided much insight into the physics of dual 
superconductivity[21-31]. Bali et al.[27, 28] presented 
the first large lattice simulation of the dual Abrikosov 
vortex adding credibility to this approach to dual su- 
perconductivity. They employed techniques to select the 
best Gribov copy [37]. In addition to improving the gauge 
fixing, this technique reduces noise coming from the ran- 
domness of false maxima. Gubarev et al. [29] found im- 
provement over the fit of Bali et al. by using the numer- 
ical solution of a lattice version of the dual GLH model. 
They concluded that dual superconductivity type lay on 
the borderline of type I and type II. Koma et al.[30, 31] 
went further and checked that the results obey scaling. 
They also argue for a decomposition of the vortex into a 
photon part and monopole part. Whereas previous work 
used a two dimensional model of the vortex, Koma et al. 
fit to a numerical solution of a three dimensional lattice 
model. Both Bali et al. and Koma et al. concluded the 
theory is described by weakly type I dual superconduc- 
tivity. 

We present results showing that the strengths of the 
dual GLH model are accompanied by weaknesses when 
applied to the interpretation of the simulation of this the- 
ory in the standard simulation window. Namely, there 
are identifiable parameters in the simulation taking dif- 
ferent values that are equal by default in the dual GLH 
model. This will inevitably increase systematic errors. 
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Also, in addition to the well-known area effect of con- 
fining magnetic currents, we show there is a significant 
perimeter effect of circulating electric currents which is 
absent in the model. Although these electric current ef- 
fects are expected to vanish for Wilson loop sources of 
very large spacial and temporal sizes, we show they can 
have an impact on the determination of the type I vs. 
type II dual superconductivity for fixed loop sizes. The 
added parameter and the electric currents are intricately 
related. 

In order to explore these issues we present an analysis 
that does not depend on the details of the dual GLH 
model. 

We regard the question of type I, type II or borderline 
to be the most important issue in these studies and we 
believe it can not be answered using analysis based on 
the dual GLH model. Perhaps the best outcome would 
be if Gubarev et al.[29] are correct that the system is 
on the borderline between type I and type II. We find 
it likely that type II behavior goes away as one goes to 
larger quark separations and better suppression of ex- 
cited states. But this is by no means definitive. Type 
I brings with it special problems which are difficult to 
address as we will explain. 

Since the primary goal is to describe the total elec- 
tric flux and the profile of the vortex as well as possible 
we consider it a top priority that our analysis satisfies 
Maxwell equations exactly for lattice averages. This re- 
moves the arbitrariness in the choice of lattice operators. 

With this background, there are in our opinion a num- 
ber of reasons to revisit the dual vortex problem. Let us 
sketch here the main issues presented in this paper: 

(i) London relation in superconductivity 

It is the lure of a general principle that draws us to 
further studies of the dual Abrikosov vortex in this paper. 
Let us first recall how this works in superconductivity [38] 
since this is the guide to the analysis of the simulation 
data. Consider an infinitely long static Abrikosov vortex. 
Recall 

J< e > =e3^*(V -ieA)</>, 

where <f) is the complex superconductivity order parame- 
ter. Further consider SGSB of U(l) via 

4> = ve^^; v constant. (1) 

Then 

= -e 2 v 2 (A- Vuj) 

and the curl of this gives the London relation 

B = -A 2 curlJ (e) , X = l/ev. (2) 

Equation (1) is the case for an extreme type II supercon- 
ductor but it is also valid deep inside a type II super- 
conductor far from a boundary where the complex phase 
degree of freedom of the order parameter dominates over 
the modulus. (The order parameter is governed by a 



'Mexican hat' potential inside the superconductor and 
by a symmetric paraboloid potential in the normal ma- 
terial. The spacial 'kinetic energy' term excites the radial 
degree of freedom in the transition region.) 

To construct a vortex with one unit of quantized mag- 
netic flux assume that the phase of the order parameter 
increases by 2tt as the azimuthal angle makes one com- 
plete path around the vortex at large transverse distances 
where the current is exponentially small. Then 

j>A-dt = <j)\7u-d£, (3) 

$ m = J B-dS = ^. (4) 

A source consisting of a hypothetical Dirac monopolc- 
antimonopole pair would generate an Abrikosov vortex 
between them containing one unit of quantized mag- 
netic flux, <f> m , which is equal to the monopolc mag- 
netic charge. All the magnetic flux from the magnetic 
monopoles is accounted for in the vortex. 

To complete the picture using Ampere's law, Eq.(2) 
becomes 

B = \ 2 {{dl + dl)B-V(V-B)}. (5) 

There are two important aspects of this that we wish 
to note. First, since V • B = p m = 0, the square root 
of the proportionality factor in the London relation is 
clearly the penetration depth in the superconducting ma- 
terial. The Ginzburg-Landau model describes real ma- 
terials hence has vanishing magnetic charge. In the four 
dimensional version, the GLH model, this is also true. 

Second is that the exterior solution to Eq. (5) is the 
Bessel function K (r/X). This solution is applicable as 
long as the superconducting order parameter ramps up 
from the normal to the superconducting asymptotic value 
on a shorter length scale than the penetration depth, i.e. 
type II. Then sufficiently deep in the superconducting 
medium the tail of the penetrating flux is described by 
Eq.(5). 

However if the penetration depth is shorter than the 
characteristic ramp-up length, i.e. type I, then there is 
no regime for which the Ko(r/\) solution is applicable. 
There is not an easy fix. The complication here is that 
the surface energy between the superconducting and nor- 
mal material flips sign[39] and an instability develops. It 
is energetically favorable to make surfaces and the behav- 
ior of the penetrating flux changes dramatically [40] See 
Tinkham[38] for a thorough treatment of the very differ- 
ent physics of the type I 'intermediate state' and the type 
II 'mixed state.' The solutions of the GL model for type I 
superconductivity in which magnetic flux penetrates ma- 
terial involves fractals on many length scales [41]. 
(ii) Exact electric Maxwell equations 
The above picture translates nicely to a U(l) lattice 
gauge theory in which a charge-anticharge pair arising 
from an Abelian Wilson loop produces a dual Abrikosov 
vortex. More specifically, Zach, Faber, Kainz and Skala 



3 



(ZFKS)[42], using a lattice Ward-Takahashi identity, 
found a specific lattice operator for electric flux giving 
the exact electric Maxwell equation for lattice averages 
independent of lattice spacing a. For example in this 
formulation the divergence of the electric field is zero ev- 
erywhere except on the Wilson line sources where it is 
equal to the charge e (see Sec. II B ). Therefore using 
the ZFKS definition the total flux in the vortex, $ e , is 
determined exactly by the charge on the Wilson loop in 
the lattice system which is of course the gauge coupling 
e. 

In 1998, DiCccio, Hart and Haymaker (DHH)[43], 
generalized the ZFKS result to SU(2) in the maximal 
Abelian gauge. This case is complicated by the existence 
of dynamical electric currents in addition to the static 
source. The extra contributions come from the doubly 
charged vector matter fields, from the effects of gauge 
fixing and from ghosts. Nevertheless by using the DHH 
definition, the total flux, &su(2) m the vortex is again an 
exact reflection of all contributions to the charge distri- 
bution. 

We consider three definitions of flux in this paper, 
ZFKS, DHH and DeGrand-Touissant[44] (DT) and find 
that to a good approximation they give the same vor- 
tex profiles, but differ by a constant scale factor. Bali et 
al.[27] did not state their definition but since they used 
the DT construction for the magnetic monopolcs, it is 
reasonable to assume they used the same definition for 
the electric flux profile. Koma et al.[31] used the DT 
definition throughout. For our DHH definition, and for 
(3 = 2.5115 we find the scale factor differs by 40% from 
the DT definition. 

(iii) Consistent magnetic Maxwell equations 

One can then also get a consistent exact magnetic 
Maxwell equation by adopting the same definition of 
flux in defining the magnetic current [45]. This is not 
the conventional procedure. The conventional one is the 
DT construction which identifies cubes with quantized 
monopoles. By using the DHH flux instead, the current 
is conserved but is not quantized in cubes. Rather it is 
smeared out among neighboring cubes. 

The magnetic Maxwell equations determine the 
solenoidal persistent magnetic currents that determine 
the shape of the vortex. 

These deviations from the more standard definitions 
are clearly tailored for this vortex problem. The lattice 
approach allows an infinite variety of definitions as long 
as they approach the same continuum limit. We argue 
in Section II how the DHH formulation might approach 
the standard discrete monopole picture in the continuum 
limit. However the DT definition is perhaps the only 
one that can describe percolating clusters[14-17] since 
their definition depends on the discretization of magnetic 
charge. We use the DT definition here for the dual vortex 
problem as a digression from the main body of this work. 
We see that it gives support to the truncation that retains 
only the percolating cluster. 

The three definitions of flux considered in this paper 



are applied to the magnetic current and again to a good 
approximation give the same vortex profiles but differ 
from each other by a constant scale factor. 

(iv) Model independent analysis 

We chose to analyze the data in such a way as to test 
the SGSB of dual U(l) independent of specifics of the 
dual GLH model. The model commits to a quartic inter- 
action. And it has vanishing dynamical electric currents. 
For gauge fixed SU (2) we determine a value of which 
satisfies a London equation asymptotically in the tail of 
the profile giving direct evidence for SGSB of dual U(l). 
That is, considering the vortex aligned along the z axis 

E z = -A^curlJ^),. (6) 

This argument does not depend on specifics of the solu- 
tions of the non-linear equations of the GLH model. 
Next consider the profile of the fluxoid[46] £ z 

£ z = E z +A 2 d (cuAJ^) z . (7) 

The leading decaying exponential behaviors E z , curl«7 m ~ 
e~ r / Xd , should cancel giving a sub-leading behavior E z + 
A^curlJ m ~ e~ r ^ d . In an extreme type II dual super- 
conductor the fluxoid would vanish exactly except on the 
axis of the vortex where the system is 'normal', i.e the 
order parameter is zero. In our case the order parame- 
ter ramps up over a finite distance and we see a second 
exponential behavior after cancelation which defines the 
length scale corresponding to the violation of the London 
relation, i.e. the dual coherence length ^[47] [48]. In this 
way we can estimate the three parameters A^, Xd and £<j 
as long as nd > 1, i.e. type II, where [49] 

Kd = —■ 

This analysis of estimating exponential behaviors is ap- 
plicable only if Xd > £d since only then can one associate 
the two exponential behaviors with the parameters Xd 
and £d properly. In other words we must have type II 
dual superconductivity for these studies to be meaning- 
ful. For type I, the identification of the penetration depth 
from the behavior Kq(t/X) breaks down because of the 
'intermediate state' which is a mixture of normal and su- 
perconducting material. Further the concept of a fluxoid 
as a sub-leading behavior breaks down. 

It is our opinion that this limitation to type II is not 
unique to our approach for reasons mentioned in subsec- 
tion (i) above. If one models the system with the dual 
GLH model and concludes type I behavior, one needs to 
show that the solutions used are capable of producing 
the serpentine surfaces[40] of the complex 'intermediate 
state'. Further one needs to address how the lattice quan- 
tum averages relate to chaotic classical solutions. 

(v) Electric currents and dual superconductiv- 
ity 

In the dual GLH model electric currents vanish every- 
where except at the location of the Wilson loop static 
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source. (Recall that the GLH model has vanishing mag- 
netic currents.) The dual of Eq.(5) implies then A<j = A^. 
In the simulation we find that there are non-vanishing 
electric currents and further as one would then expect 
\d 7^ Ad- This argues that our analysis can not be based 
on the details of the dual GLH model if we wish to learn 
more about these effects. 

We are reporting here that the (z, t) curl of the electric 
current is non- vanishing on the same plaquette where we 
measure the z component of the electric field which is 
the same plaquette where we measure the (x, y) curl of 
the magnetic current, i.e. (cmlJ^ m ^) z , out to rather large 
transverse distances. A very localized cloud around the 
sources was well-known e.g. [27, 43, 51]. But our results 
mean that this perimeter effect is much more pervasive 
than expected. For sufficiently large R and T dimensions 
of the Wilson loop, the currents and fields approach con- 
stants in z and t in the transverse midplane, resulting in 
the vanishing of the (z, t) curl of the electric current. But 
these currents are relevant in the simulation window and 
need to be understood. Determination of dual supercon- 
ductivity parameters can only be definitive if the data is 
interpreted with the correct model. 

Even if the electric currents were to vanish, we point 
out that Ad depends on the definition of flux relative to 
the definition of magnetic current through the London re- 
lation, whereas A<j is obtained from a characteristic decay 
length independent of normalization. In the dual GLH 
model these are of course equal and hence failure to be 
consistent in the two definitions will lead to a compro- 
mised fit. Hence we argue for particular definitions based 
on the physical principle of satisfying Maxwell equations. 

In Section II we review and compare three definitions 
of flux that have appeared in the literature. We argue 
in favor of particular forms of the flux and current which 
will satisfy the Maxwell equations for lattice averages. 
We show how the distinction of the two parameters A^ 
and Xd are related to the existence of electric currents 
in the vortex. The generalization of the dual of Eq.(5) 
to the four dimensional lattice problem is the key. This 
relation is exact for lattice averages only with our choice 
of definitions of flux, i.e. ZFKS for U(l) and DHH for 
SU(2), respectively and corresponding electric and mag- 
netic currents which respect Maxwell's equations. 

The numerical results are given in Sec. III. We com- 
pare the three definitions of flux. Then using our pre- 
ferred definition of flux, we measure the total flux in the 
vortex. We then implement the model independent anal- 
ysis described in sub-section (iv). We find type II behav- 
ior for Wilson loop sizes where others [27, 31] find type I 
behavior. Although these differences may be the result of 
differences in our approaches, or in the details, we regard 
the dual GLH model as the central culprit in the case of 
Koma et al. which we discuss in Sec. Ill C. 

A comparison of our work with Koma et al.[31] and 
Bali et al. [27] shows that our determination of A^ is 
essentially unaffected by the Gribov copy problem at 
(3 = 2.5115 which we discuss in Sec. Ill C. The other 



parameters could have a small dependence as reflected 
by the small Gribov copy effect reported in the tails of 
the profiles[31]. Bali et al. makes use of the London rela- 
tion as a first estimate of the penetration depth as we do 
which allows an interesting comparison possible between 
the results of these two groups and our results. 

We give examples of profiles of the electric currents 
that are contrary to the dual GLH model. 

Also as a digression from the main body of this paper, 
we show the effect of using DT definition of magnetic 
current compared to a truncated DT definition in which 
only the single dominant percolating cluster is included. 

Section IV gives our summary and conclusions. 



II. THREE DEFINITIONS OF FLUX 

Let us consider three definitions of field strength or 
flux, all agreeing to lowest order in the lattice spacing 
a. If we require that both electric and magnetic Maxwell 
equations for lattice averages of flux and current be satis- 
fied, then the specific form of the action implies a unique 
definition of flux. Using any of the alternative defini- 
tions introduces contributions, non-leading in a, that vi- 
olate Maxwell's equations. Consistency would then be 
restored only by going to the continuum limit. 

As we will see in Sec. II D below it is trivial to get an 
exact magnetic Maxwell equation as an operator relation. 
This is not possible for the electric Maxwell equations. If 
one had an exact electric Maxwell equation on a partic- 
ular configuration, updating one link can change a single 
plaquette in the electric divergence violating the equality. 
Whereas in the magnetic divergence, changing one link 
always changes two plaquettes in a compensating fashion. 
Attaining exact electric Maxwell equations is not trivial: 
the flux definition depends on the form of the action and 
the equations hold only for lattice averages. 



A. DeGrand Toussaint[44] (DT): F$ 

The first definition is that used by DeGrand and Tous- 
saint to define discrete monopoles in the U(l) theory: 

F«(n) ee ^(n)-27rn^(n), (8) 
8^(n) ee 0„(n) -e^n + iy) - 6 v {n) + 6 v {n + n), 

where # M refers to the U (1) link angle in the domain — n < 
Ofi < +7r. The integers n M „ are determined by requiring 
that — 7r < F$ < +7r. That is F$ is a periodic function 
of 6^, with period 2-7r[50]. We also refer to F$ as the 
'sawtooth' flux as shown in Fig.l. 
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plaquette angle 



FIG. 1: F$ (sawtooth) and F^V (sine) as a function of the 
plaquette angle 6* M „ 



B. Zach-Faber-Kainz-Skala[42](ZFKS): F$ 



A second definition 

F$(n) = sin^(n), 



(9) 



has the property of giving the exact electric Maxwell 
equation for lattice averages for the case of the U(l) 
gauge theory with Wilson action. Consider 

Z w (e li (m)) = sin 9 W exp ((3S) , 



S = ]T (cos^(n)-l), 

n, f-i>v 



(3=-. (10) 



The subscript of Z w {e^{m)) refers to the incorporation 
of the source into the partition function and the argu- 
ment is a variable defined as the shift of one particular 
link, 9^(m) — > M (m) + e^(m). This translation can be 
transformed away since the measure is invariant under 
such an operation. 
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w 



5e 
where 



= 



HK{f&) w -(%>) w , (ID 



(...) 



(• • • s'mdw) 



w 



(cos 6w) 

and where difference operators are defined 

A+*(n) ee *(n + e M )-*(n), 
A"*(n) = *(n)-*(n-e M ), 



where is a unit vector. 

The current {j^^ is that carried by the Wilson 

loop, normalized to take values —1,0,1, i.e. no factors 
of e. (The expectation values are superfluous for this 
observable since the shifts produce these values directly 
from the Wilson loop factor in the partition function. 
The LHS of Eq. (11) builds up these integer value from 
the lattice average.) 



C. DiCecio-Hart-Haymaker[43](DHH): F$ 

The third definition is applicable specifically to 
Abelian projected SU{2). We restrict our attention to 
the maximal Abelian gauge defined as a local maximum 
of 

R = ^tr{<7 3 ^(n)a 3 ^(n)}, 



the 



set of gauge transformations 
U — > U 9 . Taking U to be 
the stationary value, the stationary condition is given by 



{g(m) = e i«i 



Fjn[U] = 



dR[U 9 } 



daj(n) 



0. 



a=0 



The second derivatives entering in the Jacobian are given 

by 



M. 



jn;im 



(U) = 



d 2 R[U 9 



daj(n)dai(m) 



a=0 



The partition function is 

Z 9 w f -{el(m)) = J[dU]^Tr[ia 3 U w (n)}x 

e X p((3S)Y[5(F jn [U}) A FP , (12) 

jn 

where the Faddeev-Popov Jacobian is 

A FP = det\M jn , im (U)\. 

An infinitesimal shift in this partition function has the 
added complication that it violates the gauge condition. 
This can be corrected by an infinitesimal accompanying 
gauge transformation. Thus the shift in one link affects 
all links. However experience has shown that the effect 
drops off rapidly with distance from the shifted link. 
Invariance under the shift leads to the definition of flux 

sin0„„(n), (13) 

where the link variables are parameterized by 0^(n), 
(/> M (n) and 7 ([t (n) 
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and where 

Cn{n) = cos0 M (n), 
S^n) = sin^(n). 

The decomposition of an SU(2) link to U(l) gives a 
U(l) gauge link, e l6 ^, and a doubly charged matter field, 
sin^e^ 7 " -6 ") transforming at the sites[10, 11]. In the 
maximal Abelian gauge they also satisfy the vector field 
auxiliary condition making them bona fide charged vector 
matter fields as reviewed in Ref.[51]. 

This form has the property of giving the exact elec- 
tric Maxwell equation for lattice averages for this case of 
SU(2) in the maximal Abelian gauge with Wilson action 



0A- 



p(3) 
pv 



7^e)total\ 

J n ] 



(15) 



monopole with equal flux out of the six faces of the cube 
(and a Dirac string extending out from any face). Then 
the ratio of the flux out of this cube compared to 
the F^> flux gives 



6sin(27r/6) 
6(2tt/6) 



0.83. 



(17) 



On a large surface the total flux is the same for the 
two definitions. Since charge is conserved, the bal- 
ance is made up by magnetic charge in the neighboring 
cubes. We interpret this to mean that with magnetic cur- 
rents defined with F^ 2 \ the discrete monopoles become 
smeared but maintain the same total magnetic charge. 



where J, 



<f(e) total 



gets contributions from the Abelian Wil- 



son loop, the charged matter fields, gauge fixing and 
ghosts. The normalization of the Wilson loop contri- 
bution to the current in this expression is analogous to 
the ZFKS case i.e. no introduced factors of e or a. See 
Ref.[43] for details. 



D. Consistency with the magnetic Maxwell 
equation 



For the second and third cases we have a unique flux 
F^u , for i = 2,3, by requiring an exact lattice electric 
Maxwell equation. Given this definition of flux the mag- 
netic Maxwell equation is 



1 



A+ - - f( m ) 



2,3. 



(16) 



which gives a unique definition of the magnetic current. 
However the monopole current is usually taken from the 
DT definition 



1 



A+p(i) 



(This current is normalized to give monopoles with a flux 
of 27m where n is integer.) Hence if we use the conven- 
tional F^ to define the monopole current, and F^ or 
F^ respectively for U(l) and SU(2) theories to get an 
exact expression for flux in the confining string, then the 
magnetic Maxwell equation is violated. 

The electric Maxwell equation determines the to- 
tal electric flux in the confining string and the mag- 
netic Maxwell equation determines the transverse profile 
through the solenoidal currents. The only way for the cal- 
culation to be consistent with both Maxwell equations is 
to relax the usual procedure using the DT monopole def- 
inition and instead use F^ or F^ when defining mag- 
netic currents for the U(l) and SU(2) cases, respectively. 

A simple configuration will help illustrate the differ- 
ence between F^ and F^ 2 \ Consider a single DT 



E. Comparison 



Figure 1 shows a comparison between the first two def- 
initions. We plot F$) as a function of 6^„, giving a 
"sawtooth" shape. Monopoles occur as a consequence of 
Q^v crossing the sawtooth edge, giving a mismatch of 2tt 

in the flux out of a cube. The sine function, Fj$ , has 
no such discreteness and so the notion of discrete Dirac 
strings and Dirac monopoles is absent. 

However as one approaches the continuum limit the 
two definitions merge. To see this we note that with in- 
creasing (5, the Boltzmann factor supresses the plaquette 
angle to be in a more restricted effective domain in the 
neighborhood of mod 2tt. But we see from Fig. 1 that 
the sin function and sawtooth function become the same 
in this domain up to corrections of the order of a 2 . We 
expect both forms to give the standard Dirac picture in 
the continuum limit. 



The F$ definition is a modification of F^ 
volving factors of C M , the cosine of the matter fields 
in the Abelian projection of the SU(2) variables. 
Poulis[52]pointed out that has small fluctuations 
and can be taken approximately as a constant, i.e. 
(Cfj, x O) « (C,j,) x (O), where O is a generic operator. 
Hence to this approximation, we find 



(2) 



j?( 3 )w 



(Cp) 4 



(The first factor on the right cancels in the physical nor- 
malization of the flux as we show in the next subsection.) 
We will show in Sec. Ill A that for = 2.5115 this ap- 
proximation has violations of the order of 10%. 

Since C M w 1 + 0(a 2 ) in the Maximal Abelian gauge 
the fluctuations of this factor are suppressed by an extra 
power of a. As an indication of the relative fluctuations 
in Cju, our measurements of field strength have errors at 
best of the order of 1.0% whereas the error in (C^) is 
0.01%. 
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F. Physical normalization of the field strengths 

The physical dimension of the fields is obtained from 
the leading order of a small a limit. The definition of 
electric charge in the U(l) theory is straightforward and 
well-known. The Wilson loop carries the charge of the 
gauge coupling constant. 

The Abelian projected SU(2) case is not as clear cut. 
Measurement of electric flux at the site of the Wilson 
loop includes the bare charge. In addition there are con- 
tributions due to the dynamical charged matter fields, 
the effects of gauge fixing and ghosts [43]. Further an 
electrically neutral combination of the matter fields can 
not be completely disentangled from the gauge fields. We 
propose a normalization that involves the effects of these 
matter fields. 



1. (7(1) theory with Wilson Action 

Given the electric Maxwell equation, Eq.(ll) and the 
magnetic Maxwell equation, Eq.(16) (for i = 2 and with 
the Wilson action, Eq.(lO)) one arrives in the usual way 
that 



1 



and e m — 



2tt 



current jjf^ is nor- 



As mentioned in Subsection B the 

malizcd to unity on the Wilson loop. Recall that J^" 
is normalized to take the value 2ir for a DT monopole. 
The argument in subsection D above shows that the same 
normalization holds for a smeared monopole constructed 
from ZFKS flux. The normalizations are 



17(2) - 


ea 2 


(18) 


/(e) = 


— JW 

a 3 " ' 


(19) 


j(m) _ 


2tt / 1 j (m) \ 
e \ 2va? " J ' 


(20) 



Using the notation of Sec. II C we write the Wilson 
action 

/3^Tr(lV) = P^CvWCvtm + riCrim + v) x 
p p 

C v (m) {cos 6> M „(m) — 1} H 

For the maximal Abelian gauge we will make use of the 
fact that the fluctuations of are suppressed. If we take 
C M to be constant and compare with the U(l) action, 
Eq.(10), we can introduce the U(l) charge through the 
relation. 



-2 = /?(C.) 4 - 



(22) 



Then we can write 



J fll> 



j{e) Abelian Wilson loop\ 



,ea 2 {Crf 



K, 



(e) 



In the simulation for (i = 2.5115 we found (C M ) = 
0.94784(4) which gives 



e = 0.7024(1) 



(23) 



This is the charge carried by the bare Abelian Wilson 
loop. Equation (23) is only an estimate of the value of 
e because has small fluctuations and can not be con- 
sidered strictly a constant. 

Considering finally the magnetic equation, Eq.(16) for 
i = 3 we arrive at a normalization analogous to that for 
the U{1) case 



F (3) = I 

_ _l?(e) 
i 

r(m) 
"V 



m 3 



_P(3) 
4 /jz-' ' 



J 1 " 
4 *> 



(24) 
(25) 
(26) 



2. SU (2) theory with Wilson Action in the maximal 
Abelian gauge 



G. Effects due to electric currents using DHH 
definitions 



Generalizing to this case we have 

^j?( 3 )^ = ^ j( e ) Abelian Wilson loop ^ 

+ (^ e) ), (21) 



where is the sum of the three dynamical terms in 

, 7(e) Abelian Wilson loop . ,. , , 

the current and is normalized to 

unity on the Wilson loop. The expectation values are 

understood to be taken in the background of an Abelian 

Wilson loop source with gauge fixed configurations. 



Consider the classical magnetic Maxwell equations 



r (m) 



Take the curl of this 



A — t(™») _ 



1 1 , 



(27) 



= --{A-A+F Qp -A-A+F Q 
+A-A+F pff }. 



(28) 
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Using the classical electric Maxwell equation, 
1 



-A~F — /( e) 
a a P a ~ p ' 



(29) 



we get the following relation 

e apafj A-J^ = (A+ J( £ ) - A+JW) + A-iA+^ pff . (30) 

Now using the DHH form in the two currents and the 
flux, i.e. , this expression must hold for expectation 
values 

e apal3 A~ (jjT^ = 

[At (J^) A+ ( J^)) + A-U+ (Fj&) (31) 



Applying this to a vortex oriented along the z axis and 
choosing a = 3 and p = 4 

I (a,- -a 2 -(,«)) 

-l(A r Aj + A 2 -A 2 +)(F< 3 3 >) = 

-I (AJ (/<->) -AJ<J5"» 

+ 1 (A s -A+ + A ( -A+) (Ff ). (32) 

The electric current can survive in the dual GLH model 
but only as a lattice artifact and it vanishes in the con- 
tinuum limit. Further the second derivative terms on the 
RHS of Eq.(32) are designed to be as small as possible 
by the choice of the source. Assume first that the RHS 
vanishes. Then the London relation becomes 



+ A3 



(^curlj ( 



m ) 



(E z ) - \ 2 d W\ (E z ) = 0. 



This clearly identifies A d = \ 2 d as the penetration depth 
in the dual superconductor. 

However the RHS does not vanish for lattice averages 
in our simulation as we show in Sec. HI. In the standard 
simulation window, the terms on the RHS can be com- 
parable to the terms on the LHS. Hence the value A^ as 
measured by the London relation in the tail of the pro- 
file does not control the rate of transverse fall-off of the 
profile. In a dual GLH model it does. For this reason we 
choose not to rely on a fit to the dual GLH model but 
concentrate instead on verifying the model-independent 
SGSB, and estimating the three parameters A d , Xd, and 



III. SIMULATION 

Our measurements were on 208 gauge-fixed configura- 
tions on a 32 4 lattice, with (3 — 2.5115. Each update 
consisted of a 10 hit metropolis sweep and an overrelax- 
ation sweep. We made 13 runs on 16 independent nodes. 



D ropping 2000 thermalization updates (on each node), 
we made measurements on every 100th update. 

We gauge-fixed to the Maximal Abelian Gauge (MAG) 
using overrelaxation with the criterion of the average of 
the absolute value of the off-diagonal matrix clement of 
the MAG adjoint operator < 10~ 6 . 

We measured Wilson loops in which spacial links were 
fattened through 100 iterative steps by adding spacial 
staples of weight equal to the original link. In a spot 
check of the data sets, increasing the number of iterations 
from 100 to 200 and/or decreasing the weighting factor 
of the staple contribution by one half had no appreciable 
effect. Although the fattening appears to have saturated, 
the data shows a residual T/a dependence, where T/a is 
the time extent of the Wilson loop in lattice units. 

All vortex profile graphs presented here are trans- 
verse slices through the mid-plane, the (x, y) plane, on 
the quark-antiquark axis, the z axis. Fitted parame- 
ter values are calculated for 6 quark-antiquark separa- 
tions R/a = 3, 5, . . . , 13 and for 7 time separation of 
T/a = 3, 4, . . . , 9. We used Minuit for the fits and quote 
Minuit errors in the parameters determined by Ay/ 2 = 1 • 

We used a 2 F^f) throughout for the definition of flux 
in order to give the correct electric Maxwell equations. 
Similarly we used the same definition in constructing the 

magnetic current a 3 jj/™^ in order to get the correct mag- 
netic Maxwell equation. (Exceptions of course include 
graphs comparing different definitions and truncated vs. 
complete monopole loops.) 

Noise in these data is a problem. We reduced this by 
taking an azimuthal angular average over an annular re- 
gion of width = 1 in lattice units in the transverse, (x, y) 
plane, weighting the data by the fractional area overlap 
of the data point plaquette to the annulus of radius r 



<->« 



2irr J r<r , 



r<r'<(r+i) 



{■■■)da'. 



(33) 



Using this together with Stokes theorem we have a con- 
venient numerical evaluation of 

<f J-dr' = 27rr(J (r)) , 

J r ' — r 

= / (curlJ) • da* ', 

J r' <r 

= - ( (curlJ) • da', (34) 

J r<r' 

where we used the identity J (curlJ) • da' = 0. 

Since the current is a first derivative of the flux and the 
curl of the current is a second derivative, errors become 
more difficult for the latter. However by integrating an 
equation over a transverse area involving the curl we are 
back to first derivatives which are more manageable. 

Except where a fit to the data is noted, the lines in the 
graphs connect the data points. All axes on all graphs 
are dimcnsionless. For physical normalizations, one can 
take the standard value a = 0.086fm = OAAGeV' 1 . 
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gauge action summand 


gauge coupling 


quantized vortex flux unit 


GLH model 


(3 cos 6> M „ 


e = /r^ 


$m = e m = 27r/e 


dual GLH model 


/3 (d ) cos0$ 


e - 3~ 1/2 


$ e = e = 27r/e m 


U(l) theory 


cos O^v 


e u(i) - P um 


*t/(l) = e C/(l) 


51/(2) MAG theory 


~ /?S(/(2) <C M } 4 cos Oy.u H 


esu(2) = Psu{2) ( C ») 2 


< I > S(7(2) = e S !7(2) + (anti) screening 




<C M ) = 0.94784(4) 






Psu(2) = 2.5115 


esu(2) = 0.7024(1) 


$ sc/(2) = 1.72(3) 



TABLE I: Relationship between gauge coupling and quantized flux in the vortex. 
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FIG. 2: Profiles for three definitions of the electric flux 

a 2 -p34 through the transverse plane. R/a = 7, T/a = 3, FIG. 3: Same as Fig. 2 for the z component of a 4 curlj' r 
with azimuthal angular averaging. 



From this point on, we suppress the expectation value 
bracket notation, leaving it understood that we are tak- 
ing expectation values with a smeared Wilson loop. All 
flux, and currents are constructed using the DHH form 
F$ except in Sec. Ill A and III E. 
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A. Comparing flux definitions 

Recall the three definitions of flux: F^J ' . Eq. 
13). Figure 2 shows the E z profile for the three definitions 
of flux. Surveying a large variety of quark separations 
R = 3, 5,..., 13 and time extents T = 3,4, ...,9 the 
three definitions appear to differ only in the scale. We 
did not notice any significant difference in shape. The 
scale factors over this range are approximately constant 
for (3 = 2.5115. 



1 if( 3 ) 



(35) 



We note that if the small fluctuations were in fact absent 



in C M , Eq.(13), then F t 



(3) 



F 



since the C M factors 



cancel out in the normalization, Eq.(24). 

Figure 3 shows the z component of curl j( m ) con- 
structed from the three definitions of flux. The same 
observations hold here as for the flux. Note that the con- 
necting lines cross close to zero, further indicating just a 
scale factor. They scale with the same factors as the flux 
to about 1%. 

Figures 2 and 3 show qualitatively the signal for a dual 
Abrikosov vortex no matter which definition is chosen. 
For large transverse distance, e.g. in this case r/a > 5, 
the London relation holds, i.e. the tails can be arranged 
to cancel, 



fluxoid : £ z = E z + A? ( curlJ (r 



0. 



(36) 



The integral of the (curU^ m -') 2 over the whole plane van- 
ishes hence the integrand must change sign. For an ex- 
treme type II dual superconductor, the profiles match 
and cancel everywhere except for a delta function con- 
tribution at r = 0, giving the sign change needed for 
vanishing of the surface integral. In our case the scale 
defining the breakdown of the London relation is the dual 
coherence length £<j. 
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Figures 2 and 3 also show that by varying the rela- 
tive definition of the two quantities one would change 
the value of A^ without changing the value of Xd leading 
to systematic errors in the analysis in a model in which 
the two parameters are forced to take on the same value. 



B. Total electric flux 

We first consider the total flux integrated over the 
transverse mid plane since these data generally have the 
smallest errors over the range of the 42 data sets. 

In Fig. 4 we plot the integrated flux and fluxoid in the 
interior of the disk r' < r 



E z da'; 



r'<r 



£ z da'. 



The integrated fluxoid reaches the asymptotic value at 
smaller r/a than the integrated flux because it corre- 
sponds to the ramping up of the order parameter, not 
the penetration of the flux. They must have the same 
asymptote because the surface integral of the fluxoid over 
the whole transverse plane vanishes. 



SU(2) 



j EZ ^ = l 



E, da' . 



Figure 5 shows the total flux obtained from the radial 
asymptote of the fluxoid as a function of R/a and T/a. 
This shows that for each R/a the data goes to an asymp- 
tote for large T/a, at least for those values of R/a where 
the errors do not mask the behavior. With this as a guide 
we do an exponential fit to extract &su(2) f° r infinite T/a 
values, based on one term in the transfer matrix. 



5(7(2) 



= A + Be 



-CT/a 




FIG. 

R/a 



4: Flux and fluxoid integrated over a disk of radius r/a. 
= 7, T/a = 3. 
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FIG. 5: Total electric flux in mid transverse plane. Within 
each group labeled by R/a, the 7 ticks correspond to T/a = 
3, 4,. ..,9. 



Further we do a global fit, using all the data available 
to arrive at a value of total flux and we then test it for 
independence of R. To do this we take the above form 
and allow an extra dependence linear in R/a in the pa- 
rameters B and C giving 



®su(2) = A+(B + B'R/a)e 



-(C+C'R/a)T/a 



The results are shown in Fig. 6. We plot the value 
of A as a function of the largest R/a included in the 
fit, -Rmax/a- Similarly, the family of curves shows the 
dependence on T max /a. The value of $su(2) determined 
this way using all 42 data sets is 



SU(2) 



-I 



E z da = 1.72(3), x7 d - f - = 1-0. 



As an alternative to this extrapolation, we can also just 
take a few data points from Fig. 5 for values of R/a and 
T/a large where the value seems to be stable but small 
enough so that statistical errors are manageable: 



R/a = 7 $ 
T/a = 6 { R/a = 9 $ 
R/a = 11 $ 



SU{2) 
SU(2) 
SU(2) 



1.76(6), 
1.75(9), 
1.81(9). 



Table I summarizes the relationship between the gauge 
coupling constant and the quantized vortex flux for the 
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10 

R max/a 



FIG. 6: Extrapolation to T/a 
transverse flux based on data 
3, 5, . . . R m ax/a for 4 values of T ma: 



= oo of the integrated 
for spacial separations 
/a. 



interpretation is that lattice average describes a vortex 
that can go in either direction and that the measurement 
on one side is a superposition of and 1 units of flux. 

Therefore since a superposition is obviously possible, 
we think it is best to regard e m — 2tt/^su<2) as a tenta- 
tive result subject to a better understanding of the dual 
GLH model or a generalization of it to include the electric 
currents. 

We wish to note some interesting aspects of Fig. 5. 
For each value of R/a, $su(2) increases with T/a. The 
Wilson loop carries an electric current and the dynami- 
cal currents circulate in the same sense. This is evident 
when looking at a fixed time slice since there is antis- 
creening in the divergence of the current. Since increas- 
ing T/a suppresses excited states, we must conclude that 
the electric current associated with them circulates in the 
opposite direction than in the ground state. Further we 
note that the saturation of the anti-screening cloud as a 
function of R/a is also evident at fixed T/a. We take the 
case T/a = 7 which has manageable errors for increasing 
R/a. In that case the anti-screening already appears to 
be saturated for R/a = 5. 



C. Three parameters of the vortex 



four cases considered here. For the U(l) and SU(2) cases, 
the Wilson loop carries the charge £u(i) an d &su(2) re- 
spectively. For the (7(1) case, the quantized unit of flux 
is also &u(i) an d so the elementary charge produces ex- 
actly one unit of quantized flux in the vortex. In the 
SU(2) case however there is a dynamical charge distri- 
bution generated by the source exhibited by the large 
value of e dynamical = § S u(2) - e SU (2) = 1-02(3). (Since 
Gauss' law is satisfied exactly in this formalism, &su(2) 
measures exactly the total charge on each side of the 
transverse plane.) 

If we are correct that the proper interpretation of the 
simulation data implies that there is a dual gauge theory 
operating, fixing the quantization of flux to one unit, 
then from lines 2 and 4 in Table I we conclude that the 
fundamental unit of flux in that dual theory is &su(2) = 
1.72(3), implying that the gauge coupling constant in 
that dual theory is e m — 2n/& SU ( 2 )- 

Is the assumption of one unit of flux correct? Could it 
be that the correct interpretation of our data is a super- 
position of vortices with multiple units of flux? 

The (7(1) example suggests that this is indeed 
possible [54]. Consider a lattice with a small extent in the 
z direction, i.e. along the alignment of the vortex. Con- 
sider the case in which the distance between the sources 
is the same in both directions around the torus. And take 
a large enough time extent of the Wilson loop so that the 
spacial direction of the flux is not biased in one direction 
or the other. We will obviously find half a unit of flux 
in each direction. We are describing a system in which 
the vortex can occur in either of the two directions. Our 
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FIG. 7: Profiles of electric flux and fluxoid integrated outside 
a disk of radius r/a. Fits corresponding to Eqs. (37) and (38) 
also shown. R/a = 7, T/a = 6. 

Surveying the quality of the data for all 42 data sets 
we arrived at the interval to identify the tails. We chose 
the 10 points: n/a = 3.5, 4.0, • • • ,8.0 for x 2 fits. 

For the flux, we search for the best values of A and 



f a 2 E z da' = A^fr~fae-^l^. 

Jr'>n 



(37) 
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FIG. 9: Fitted parameters Ad/a based on Eqs. (38). Within 
each group labeled by R/a, the 7 ticks correspond to T/a = 



FIG. 8: Fitted parameters Ad/a and £d/a based on Eqs. (37) 
and (38). Within each group labeled by R/a, the 7 ticks 
correspond to T/a — 3,4, ... ,9. 



Figure 7 shows the integrated flux profile along with a 

x 2 fit. 

In our analysis, we have used exponential tail forms 
appropriate for two dimensions. The London relation in 
this case leads to exponential behavior K n ((), see e.g. 
Rcf.[38], where K v is the modified Bessel function. We 
make use of the first term in the divergent asymptotic 
expansion 

*.<o-^( 1 + ilpi + ...). 

The two dimensional choice requires roughly r/a < R/2a 
which can not be achieved in the tail of the profile. For 
example for R/a = 8 and a transverse distance r/a = 8, 
the distance of the observation point to the sources is 8.9 
in lattice units, compared to a distance of 8.0 from the 
vortex, which is clearly not a two dimensional problem. 
Nevertheless it is reasonable at intermediate transverse 
distances and becomes correct as R — > oo and for large 
transverse distances. And it allows a close comparison 
with Bali et al.[27] 

Since we use expressions for the profile functions that 
have been integrated over the area, consider 

poo 

/ C"^+l(CK = rf v K v {j]), forz/=-l and 
J-q 

tf v _i(0-#„+i(0 = —K V {Q, for^ = 0. 
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d 

FIG. 10: Scatterplot of \d/a vs. Ad/a. 



Hence 

As explained above the choice of integrating the flux re- 
duces noise in the signal. 
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Next we do a \ 2 fit of the fluxoid in a similar fashion 




B^V~fae- Ti l u . (38) 

determining A d , B, and £<j. Figure 7 also shows the sub- 
leading behavior after canceling the leading large r be- 
havior between E z and (curlJ*-" 1 ')^, i.e. the fluxoid pro- 
file. Also shown is the fit. 

This demonstrates the heart of our method. Through 
the fit we find a proportionality constant such that the 
London relation is satisfied asymptotically by the pen- 
etrating tails of the flux and the curl of magnetic cur- 
rent respectively. Further the violation of the London 
relation is characterized by a second exponential at a 
shorter length scale. Even without a fit, the data points 
alone support the cancellation inherent in the construc- 
tion. For specificity, we chose this operational defini- 
tion of £d with the understanding that it sets the scale 
but there can be a factor close to one relating it to the 
Ginzburg-Landau definition [49]. Further the borderline 
between type I and type II also has a factor close to one 
[51]- 

Using the parameters Ad and £<j specifically defined 
by Eqs.(37) and (38) then the condition Eq,(8), k = 
^d/£,d > 1 unequivocally signifies type II behavior. Un- 
der that condition, there exists a region, possibly asymp- 
totically, where the London relation is satified and that 
eliminates type I as a possibility. 

Out of the 42 x 2 cases, X 2 /d.f. < 1.0 in the fits except 
for a few cases at small R/a and small T/a. The results 
of the fits are shown in Figs. 8, 9 and 10. 

Note first in Figs. 8 and 9 that there are less than 
42 data points for each quantity. This is because we 
set an inclusion criterion Ad > l.l£d- The reason is that 
the method of isolating the two exponential forms breaks 
down as the parameters approach each other. When they 
are close, large errors develop in trying to separate the 
two behaviors. Hence the 10% cut. The absence of a data 
point indicates we tried and failed to produce a fluxoid 
which indicates that the parameters appear to be driven 
close to or into the type I regime. Therefore every value 
of R/a and T/a with a data point indicates that we are 
able to separate the leading and subleading exponential 
forms which implies unambiguously type II behavior. 

For R/a — 5 and 7 and for all values of T / a except one, 
Fig. 8 shows that we were able to construct a fluxoid with 
the pattern shown in Fig. 7. Whether one accepts our 
parameterization of the tails or not, this is strong model 
independent evidence that for these cases, the system 
behaves like type II. On the other hand for the values 
R/a = 6 and T/a — 6, fitting to the dual GLH model, 
Koma et al.[31], find type I behavior. Because of differ- 
ences in the number of smearing steps and the value of 
the smearing paramenter, we should not compare the two 
results for individual values of T/a. However by calculat- 
ing at many values of T/a and bracketing their value of 



R/a we are sidestepping that problem in a comparison. 

This conflict could probably be checked out without 
any fitting of profiles or any commitment to a partic- 
ular method of interpreting the data. Simply take the 
data points of the curl of the magnetic current profile 
and multiply by an arbitrary constant A^ and scan its 
value looking for cancellation with the data points of the 
electric profile. A log plot will reveal if the search turns 
up a subleading exponential in the difference data points. 
If so it can only mean type II, if not then it is borderline 
or type I. 

We have also repeated the whole analysis for a generic 
exponential fit to the tails in Eqs. (37,38) without regard 
to the power behavior multiplying th e ex ponential, i.e. 
without the two-dimensional factor Ur/a. The pattern 
of results for Ad and ^ were quite similar to those in 
Fig. 8 but the values of the two parameters tended to 
be larger. The present choice gave smaller errors and a 
more coherent scatterplot, Fig. 10. 



1. Comparison to Bali et al.[27] 

Our results allow a rather direct comparison with Bali 
et al. [27] . It is interesting to compare their methods and 
best estimate of A d /a to ours. They chose T/a — 6, 
R/a = 8. Here are their steps in order: 

(1) London relation in the tail of the flux, Ad/ a = 1.82(7). 

(2) Include physics of dual GLH 

(a) but fitting only the flux, A d /a = 1.84(8). 

(b) Including the magnetic current, Ad/a = 1.99(5). 

(c) An alternative to (b), Ad/a = 1.62(2). 

(3) Best estimate: A d /a = 1.84^. 

We regard step (1) as an essentially correct definition of 
the penetration depth and with minor differences their 
methods agrees with ours[53]. Our value e.g. Ad/a — 
1.96(3) is for T/a = 6, R/a = 7 . Koma et al.[31] find 
A d = 2.22(1) for T/a = 6, R/a = 6. 

Step (2a) is a refinement acknowledging that there may 
be a non-zero coherence length detectable at intermediate 
distances in the tail. 

In the succeeding steps (2b) , (2c) and (3) we differ since 
they rely on the dual GLH model. The fitted parameter 
representing the penetration depth was found to take two 
different values statistically inconsistent with the first 
determination and each other when fitting two different 
ways. A similar inconsistency would be expected based 
on the thesis of this paper if one tailored the fit to fa- 
vor the Ad definition of penetration depth vs. one that 
favors the A d definition. The quoted systematic error in 
the parameter was triple the original statistical error in- 
dicating to us trouble since the first approach was correct 
in our opinion. They attributed the systematic error to 
lattice artifacts coming from the fitting t definition. In 
our opinion we regard this as evidence for problems with 
the dual GLH model interpretation of the data. 
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We may or may not differ on type I vs. type II. Brack- 
eting their value T/a — 6, R/a — 8 we find type II for 
the same time extent and R/a = 7, and our method 
breaks down for our next value R/a = 9. They find type 
I. Given this finding, ironically it negates their starting 
point, step (1), in using Ko(r/Xd) since that form is not 
applicable in the case of type I. (It is not clear to us if 
the succeeding steps depend on step (1).) 



2. The effect of Gribov copies 

There is more we can say. Bali et al.[27] used the meth- 
ods of Bornyakov et al. [37] to find the best maximum of 
the gauge fixing algorithm. This allows the use of suc- 
cessive estimates to look for gauge copy dependencies in 
results. On the other hand we did a simple overrelaxation 
gauge fixing algorithm. The question arises whether this 
difference of our methods with Bali et al. (step(l)) could 
lead to differing results in fitting parameters. It turns 
out that the results are insensitive to the maximum se- 
lection criteria of Rcf . [37] . To see this consider the upper 
left hand graph in Fig. 4 of Koma et al.[31]. We see that 
the flux profile is remarkably insensitive to Gribov copies. 
There is a very small Gribov copy dependence in the core 
of the vortex and it appears to entirely go away in the 
tail in the range where we and Bali et al. determined the 
tail, r/a > 3.5. If we were to take the data from Ref. [31] 
and use our model independent methods we would obvi- 
ously have to find our A<j parameter insensitive to copies 
as would Bali et al.[27] (in their steps (1) and (2a), using 
the designations in the previous subsection). 

There is another parameter we call A^, the square of 
which is the proportionality factor between the tails of 
the flux and the curl of the magnetic current. The upper 
right hand graph in the same figure, Fig. 4, Ref. [31] , shows 
there is detectable but small Gribov copy dependence for 
the magnetic current in the tail, for r/a > 3.5 where we 
work. Using their data and our analysis, that implies A^ 
could have a small Gribov copy dependence. But if their 
analysis gives a penetration depth with measurable copy 
dependence then it must come from an artifact of the 
model since two independent parameters, A<j and A<j are 
forced to be equal by the model. 

If the dual GLH model were the correct model of the 
data then the fitted parameters would have to give a 
faithful representation of the physical observables appar- 
ent in the tails. And this would be a very strong test of 
the model if the parameters truly represented the tails. 
Unfortunately there are no log graphs of profiles in the 
literature where one could independently verify the fits 
by eyeball. 

Next, consider the effect of copies on the coherence 
length. Indeed our £<j may very well have Gribov copy 
dependence, though small, for the same reason that A^ 
might show this dependence. 

That is not to say the the refinement of using the 
Rcf. [37] is not necessary. For higher values of (3, Koma 



et al.[31] have shown that the Gribov copy dependence 
is more pronounced. 



3. Other features 

Some general features of Fig. 8 are worth noting. In- 
creasing T/a drives the and £d together. With limited 
statistics we can not take R/a large enough to make a 
definitive statement as to whether type II behavior per- 
sists. 

The results for A<j are shown in Fig. 9. We also give 
a scatter plot of A^ vs. A<z, Fig. 10. There is no clear 
tendency for them to be equal as they are in the dual 
GLH model. 



fitting range, r/a 


[3.5, 8] 


[2.5, 8] 


[3.5, 10] 


no. of pts. 


10 


12 


14 


A d /a 


1.96(3) 


1.92(2) 


2.00(3) 


X 2 /d-f. 


0.40 


0.56 


0.63 


U /a 


1.21(20) 1.28(11) 


1.24(16) 


A d /a 


1.50(11) 


1.47(9) 


1.49(8) 


X 2 /d.f- 


0.03 


0.06 


0.14 


®SU(2) 




1.72(3) 




X 2 /di. 




1.0 





TABLE II: First data column: parameters determined from 
T/a — 6, R/a = 7 using standard interval. Second data col- 
umn: sensitivity to lower cut. Third data column: sensitivity 
to upper cut. The value of $su(2) does not involve cuts. 

In Table II we show the sensitivity of the parameters 
for a particular Wilson loop size to varying the fitting 
interval. Note that the second row of values of \ 2 /df for 
the fluxoid fit are exceptionally small. From Fig. 7 we 
see that the fit to the points of the fluxoid is excellent but 
the error bars appear to be too large given the smooth 
behavior of the data points. We attribute this to error 
propagation in forming the fluxoid where there might be 
correlations in the two canceling contributions. 



D. Electric currents 

Figure 11 shows an example of the effect of electric 
currents in this problem. The four graphs correspond to 
the four terms in the identity (a rewriting of Eq.(32) in 
dimensionless form and supressing the expectation value 
bracket): 

= a 3 (A-,4 m) A 2 " jM) a 2 (Ar A+ + A 2 "A+) F%> 
W (A 3 +Jf - A+ J^) - a 2 (A 3 "A+ + A 4 -A+) jrg>. (39) 

The curl of the magnetic current is in the (x, y) plane 
and the curl of the electric current is in the (z,t) plane. 
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FIG. 11: The four contributions to the identity, Eq.(39). The 
sum of the four graphs vanishes identically. R/a — 7, T/a = 
3. 
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FIG. 12: Same as Fig. 11 except T/a = 6 



FIG. 13: The (z, t) curl of the electric current evaluated on 
the vortex axis. Within each group labeled by R/a, the 7 
ticks correspond to T/a = 3, 4, . . . , 9. 



All four terms live on the same plaquette, which is the 
same as the E z plaquette. 

In a continuum dual GLH model, the third term, the 
electric current term, vanishes. Further if we consider an 
infinitely long static vortex in the dual GLH model, then 
the fourth term also vanishes. 

Both the third and fourth terms are significant and 
remain so over a wide range of i?'s and T's presented here. 
The discrepancy between A<j and Ad further support the 
presence of the third and fourth terms over this range. 

As we go to larger values of T/a the behavior of the 
electric currents is more complicated as shown in Fig. 12. 
In particular there are wiggles in the electric current term 
in the tail. There are examples of sign changes for other 
values of R/a and T/a. It is beyond the scope of this 
paper to delve into such features since we do not have a 
well motivated model of these currents. Rather we wish 
to point out here the evidence that the electric current 
persists for the range of R/a's and T/a's studied here. 

Fig. 13 shows the value of this term on the axis. For 
R/a = 3 the curl term on axis picks up the large an- 
tiscreening contribution adjacent to the source which is 
discussed elsewhere[27, 43, 51]. At R/a = 5 and larger, 
this value drops. But it persists and can change sign. 

The (z, t) laplacian term also persists. It measures how 
well we are in the regime of large R/a and T/a. The non- 
vanishing of this term contributes to the non-equality of 
Ad and A^. Such a term will arise in the three dimensional 
dual GLH model used by Koma et al.[31]. It would be 
interesting to see how well it fits this term. From our 
point of view we need a three dimensional model that 
does not commit one to the dual GLH model if we are to 
resolve the issues raised here. 

If one were to drop all expectations about dual super- 
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conductivity and simply ask in an unbiased way how well 
the dual GLH model fits the simulation data, then one 
ought to do a \ 2 fit that includes all four of these terms. 
They all live on the same plaqucttc. Then one would 
be truly testing the validity of the dual GLH model. 
Of course the result would be that the contribution to 
X 2 involving the electric current would diverge as statis- 
tics improved. And the contribution involving the (z, t) 
Laplacian would diverge if the source in the model did 
not match the source in the simulation. Measuring the 
latter term offers a good test of how well one matches the 
idealization of an infinitely long static string. 

E. Truncated monopole loops 




r/a 



FIG. 14: The effect of truncation of monopole loops, keeping 
only the percolating cluster. 

Finally we give a result that is disjoint from the 
body of this paper. We looked at monopole loops us- 
ing the standard DT definition, i.e. constructed us- 
ing F$ . In this case there is a distinction between 
connected and disconnected loops. Earlier studies have 
shown that there is a single large connected percolating 
cluster of monopole currents that dominates confinement 
physics[14-17]. There is a sharp distinction between this 
percolating cluster and the large number of very small 
loops. 

We confirm this result in our configurations and find 
that the big cluster contains ~ 40% of the current. We 
compared our calculations of the current with the trun- 
cated version, dropping ~ 60% of the contribution to 
the current coming from the small loops. Our results 
are shown in Fig. 14. This gives the integrated curl 
J r , (curl.7) ■ da' as a function of r. The truncation has 
a very minor effect in the core and essentially no effect 



in the tail which means that it has a minor effect in de- 
termining the parameters of the vortex. 



IV. SUMMARY AND CONCLUSIONS 

We regard the most interesting open question in stud- 
ies of the dual Abrikosov vortex is whether the system 
behaves as type I, type II or on the border line. We 
have cast doubt on the use of the dual GLH model as a 
reliable way to analyze simulation data. There are no dy- 
namical electric currents in the dual GLH model for the 
same reason there are no dynamical magnetic currents 
in the GLH model. Hence it is wrong, but how wrong? 
We argue that these currents are important enough to 
cast doubt on the prior efforts to establish the type of 
superconductivity. 

How reliable is our model independent approach in 
trying to answer the question? The definition of the 
GL coherence length £ and the GL borderline parame- 
ter k w A/£ are not easily divorced from the GL theory. 
And we see in real world superconductivity there are var- 
ious ways to set the coherence length scale, but they are 
related factors close to 1, sec Tinkham[38] and [49,51]. 
We have sidestepped these issues that involve the dual 
GLH model and instead defined the parameters based 
directly on the behavior of the flux and fluxoid, See Sec. 
Ill C . 

Lets put off that question and try to answer a simplicr 
question. What is the apparant type I/II behvior when 
the system is probed by a specfic Wilson loop? We have 
proposed a way to answer it. 

Consider Fig. 7. We arrived at the fluxoid profile by 
finding a value of A^ such that A^ x curl of the magnetic 
current cancels against the electric flux in the sense that 
the difference has an exponential behavior with a shorter 
characteristic length. Ignore the fits in the graph and 
the details of the functional forms. One can arrive at 
this just by stepping the value of A^ until this value is 
found. It doesn't even require determining or 

The data points alone in Fig. 7, makes it very clear 
that the London equation is satisfied asymptotically for 
large r/a. This is unequivocal type II behavior for this 
R/a = 7, T/a — 6 Wilson loop source. The London equa- 
tion is irrelevant to type I behavior because the ramp-up 
of the order parameter occurs over a larger distance than 
the penetration of the flux. 

Fig. 8 shows 25 Wilson loop sizes for which that con- 
struction was successful. There may be more but we 
made a conservative cut in the data as described in Sec. 
Ill C. There is clearly room for improvement in our pro- 
cedure to identify the fluxoid. As we go to larger loop 
sizes and better suppression of excited states, the system 
seems to be moving toward type I behavior. 

Koma et al.[31] reported type I behavior for R/a = 
6, T/a = 6. We report 13 loop sizes that bracket their 
loop size for which we are successful in the fluxoid con- 
struction. So we would surmise that had we calculated 
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R/a = 6 we would most likely report type II for all our 
values of T/a with the largest value possibly in question. 
Bali et al.[27] reported type I or R/a = 8, T/a = 6. We 
did not calculate this case or bracket it well enough to 
speculate what we would find. 

Gubarev et al. [29] reported that the system may lie on 
the boundary. This would clearly be the most satisfy- 
ing outcome theoretically since it implies there may be 
something universal happening. 

Given that the total flux in the vortex and the shape 
of the vortex are intimately connected to the electric and 
magnetic Maxwell equations, we considered it a high pri- 
ority to use a definition of flux that guarantees the equa- 
tions are satisfied for lattice averages. A bonus was an 
identity relating the magnetic and electric currents which 



was exact only for our choice of definitions of flux. This 
identity, Eq.(39) was central in elucidating the problems 
associated with the dual GLH model. 
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